Quantum imaging with incoherent photons 
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We propose a technique to obtain sub- wavelength resolution in quantum imaging with potentially 
100% contrast using incoherent light. Our method requires neither path-entangled number states nor 
multi-photon absorption. The scheme makes use of N photons spontaneously emitted by N atoms 
and registered by N detectors. It is shown that for coincident detection at particular detector 
positions a resolution oi \/N can be achieved. 
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In Young's double slit experiment (or in a Mach- 
Zehnder interferometer) the probability G^^\y) to de- 
tect a photon at position r results from the interfer- 
ence of the two possible paths a single photon can take 
to reach the detector. This is expressed by the state 
IV^(l)) = I/V2 (|l)t/|0)L + |0)c/|1)l) where the subscript 
L (U) denotes the lower (upper) arm of the interfer- 
ometer. Variation of the detector position leads to a 
modulation of the form G^^\r) ex 1 + cos(^(r), where 
6{r) = kd sin 0{r) is the optical phase difference of the 
waves emanating from the two slits and /c, d and ^(r) 
are the wavenumber, slit separation and scattering angle, 
respectively. Obviously, the fringe spacing of the modu- 
lation (in units of dsm6{T)) is determined by the optical 
wavelength A, in correspondence with the Rayleigh crite- 
rion 

Quantum entanglement is able to bypass the Rayleigh 
limit 0, S, i, i, 0, 0, 0, i, H Consider for 

example the path-entangled A^-photon state |'^(A^)) = 
1/V^{\N)u\0)l + \0)u\N)l). Because the A^-photon 
state \N) has TV times the energy of the single pho- 
ton state |1) in a given mode it accumulates phase N 
times as fast when propagating through the setup. This 
gives rise to an N photon absorption rate of the form 
G^^\r, •••,!*) oc 1 -h cosN6{r) exhibiting a fringe spac- 
ing N times narrower than that of G^^\r) [^]. This gain 
in resolution can be fruitfully applied for a wide range 
of applications, e.g., in lithography [3, Hf, microscopy [s], 
spectroscopy [9] and even magnetometry [l^ . In order to 
implement this A^-fold increase in resolution commonly 
an entangled state of the form \ip{N)) in combination 
with a non-linear medium sensitive to A/'-photon absorp- 
tion is needed [lH . 

In this letter we propose a different scheme to achieve a 
resolution of X/N involving neither of the above require- 
ments. In what follows we will apply this scheme in the 
context of microscopy. The method employs N photons 
spontaneously emitted from N atoms subsequently de- 



tected by N detectors where by means of post-selection 
it is ensured that precisely one photon is recorded at each 
of the N detectors. We demonstrate that in this case, 
for certain detector positions r2, . . . ,rAr, the A^th order 
correlation function as a function of ri takes the form 
1 + cosA^(5(ri), resulting in a phase modulation with a 
theoretical contrast of 100% and a fringe spacing deter- 
mined hy X/N. As with path-entangled number states, 
this corresponds to an A^-fold reduced fringe spacing com- 
pared to G'^^^ (r) while keeping a contrast of potentially 
100%. Hereby, only tools of linear optics are employed 
as a single photon is registered at each detector. 

To understand this outcome let us consider A^ identical 
two-level atoms excited by a single laser 7r-pulse. After 
the spontaneous emission the N photons are registered 
by N detectors at positions ri,...rAr. For the sake of 
simplicity let us consider coincident detecti on 1 121 ]. In 
that case the Nth order correlation function Il3[ can be 
written (up to an insignificant prefactor) as [ij] 

GW(ri, . . . , riv) = {D\ri) . . . D^{rM)D{rM) . . . D{r^)) 

(1) 

where 



-i/cn(ri) -R-Q, 



(2) 



Here n{ri) = r^/r^ stands for the unit vector in the di- 
rection of detector i, the sum is over all atom positions 
RcK, k = c^o/c, with coo the transition frequency, and 
o'a = \9)a{e\ is the lowering operator of atom a for the 
transition |e) \g). 

For all atoms initially prepared in the excited state |e), 
we obtain from Eqs. ([T]) and ([2]) 



GW(ri,...,r^) 



|7(ri,...,rjv)|^ 



(3) 
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FIG. 1: (Color online) Atomic arrangement and detection 
scheme: N identical two-level atoms at Ri, . . . , Rat sponta- 
neously emit N photons after excitation by a laser pulse. The 
photons are recorded in the far field by N detectors positioned 
at ri, . . . , Fat. The figure exemplifies the case = 4 (for ad- 
ditional symbols see text). 



where [15| 

7(ri,...,rAr) = 



N 



N 



]-Jg-z/cn(reJ.R«^ (4) 

ei,...,eAr = 1 «=1 

ei / . . . 7^ CAT 

Equations (|3]) and (|4]) show that ^^^^(ri, tat) results 
from the interference of N\ terms, associated with all 
possibilities to scatter N photons from N identical atoms 
which are subsequently registered by N detectors. 

To simplify further calculations let us consider the case 
of N equidistant atoms. By choosing the origin of the 
coordinate system in the center of the atomic chain, we 
have 



(5) 



with u the unit vector along the chain axis, d the inter- 
atomic spacing and ja = —{N — l)/2, . . . , {N — l)/2 for 
a = 1, . . . , A" (see Fig. 1). By defining 



^(r^) = kdn{ri) ■ u = kdsinOi 



(6) 



where Oi is the angle between n(r^) and the direction 
normal to the atomic chain (see Fig. 1), we find 



GW(ri,...,r^) = ^[^cos(j.,5)]' 



(7) 



Here, j is the vector of the distances of the atoms from 
the origin in units of d: 



S is given by: 



(ji,---,jiv), 



((^(ri),...,(5(r^)). 



(8) 



(9) 



and the sum in Eq. ([7j) is over the A^! permutations of 
the j components. 

Due to the symmetry of the configuration, the func- 
tion G*^^^ (ri, . . . , ttv) contains A^!/2 cosine terms, each 



oscillating in general with a different spatial frequency. 
Obviously, the complexity of the expression rises rapidly 
with increasing atom number A". However, if the A" de- 
tectors are placed in such a manner that all terms in 
Eq. ([7j) interfere to give a single cosine, one is left with 
a modulation oscillating at a unique spatial frequency. 
This occurs in the following case: for arbitrary even N 
and choosing the detector positions such that 



5{r2) = -S{v,) ., 



27r 



5(r3) = 5(r5) = ... = (5(r^-i) = - 
5{r,)=5{re) = ... = d{r^) = --, (10) 

the A^th order correlation function G^^^ as a function of 
detector position ri reduces to 



GW(ri)=Ayv [l + cos(iV5(ri))], 



(11) 



where A at is a constant which depends on A". For ar- 
bitrary odd A" > 1, and choosing the detector positions 
such that 



S{V2) = -5{ri) , 
6{rs) = (5(r5) = . 

,5(r4) = S{re) = . 



S{vn) 
6{rN-i) 



2tt 
N+1 



27r 
' N+1 



, (12) 



the A'th order correlation function G^^-* as a function of 
ri reduces to 

GW(ri) =An[1 + cos((iV + 1) Sir,))]. (13) 

According to Eqs. (pTj) and ([13]), we are able to ob- 
tain for any A^ a correlation signal with a modulation 
of a single cosine, displaying a contrast of 100% and 
a fringe spacing determined by X/N {X/{N + 1)) for 
even (odd) A^. We note that due to the limited detec- 
tor sizes and the dipole emission pattern of the spon- 
taneously emitted photons only a subset of all emitted 
photons will be recorded. However, in contrast to us- 
ing maximally path-entangled A'-photon states we are 
able to avoid in this scheme both the necessity to gener- 
ate a state of the form |'^(A^)) and the need to detect a 
multi-photon absorption signal pj|]. We emphasize that 
as the photons are produced by spontaneous decay the 
interference signal is generated by incoherent light. We 
stress further that a fringe contrast implied by Eq. ([TT]) 
or Eq. ([131) proves the underlying quantum nature of the 
process [ItI, lisl . [l9| . The quantum character is gener- 
ated by the measurement process after the detection of 
the first photon. In fact, just before the detection of 
the A'th photon, the atomic system is in an A^-particle 
ly-state with one excitation [20[ • The non-classical char- 
acteristics of our scheme are thus another example of 
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detection induced entanglement of init ially uncorrelated 
distant particles [11, [H, [13, [H, [26| . 

To exemplify our method, let us consider the sim- 
plest situation, i.e., the case of N = 2 atoms. With 
j = (— |, + we obtain from Eq. ([7j) 

G(2)(ri,r2) = ^ [1 +cos(5(ri) - Sir^))] . (14) 

Obviously, the modulation of the G*^^^ (ri, r2)-function 
depends on the relative position of the two detectors (see 
Fig. 2): for S{t2) = ^(1*1) the second order correlation 
function is a constant, whereas for fixed ^(r2) the two 
photon coincidence as a function of S{ri) exhibits the 
same phase modulation and fringe spacing as G^^^ (r) in 
Young's double slit experiment. However, the increased 
parameter space available for the detector positions in 
case of two detectors allows also to pick out the relative 
orientation S{r2) = —5{ri). In this case we get 



G(^)(ri) = -[l + cos(2.5(ri))], 




G(2)(ri) 
2A,i. _ 



6{r2) = 5(ri) 



G(2)(ri) 



(15) 
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47r S{ri) 
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1 ^(i"2) — const 


0- 
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exhibiting a phase modulation as a function of ri with 
half the fringe spacing of G^^^ (r) while keeping a con- 
trast of 100% (see also [HI). Note that the assumed con- 
dition for the direction of emission of the two photons, 
i.e., ^(r2) = — (5(ri), corresponds to a space-momentum 
correlation of the photons identical to the one present in 
spontaneous parametric down conversion ^, 0, [2^ . 

In the case of the fourth order correlation function 
G^^^(ri, r2, ra, r4) for four equidistant atoms, and by 
placing the detectors according to Eq. ([TO]) (see Fig. 2), 
one finds 



[l + cos(4<5(ri))] 



(16) 



Obviously, G^^) as a function of ri exhibits a modulation 
of a single cosine with a contrast of 100%, in this case 
with a fringe spacing determined by A/4. 

As an example, let us apply our scheme in the context 
of microscopy. From Abbe's theory of the microscope we 
know that an object can be resolved only if at least two 
principal maxima of the diffraction pattern are included 
in the image formation [l|. According to this criterion 
the use of the first order correlation function G^^^(ri) for 
imaging N equidistant atoms allows at best to resolve an 
interatomic spacing equal to A [l|. Indeed, if each atom 
is initially prepared in the state \(j)) = -^(Ig) + |e)), we 
get from Eqs. ([T]) and ^ 



G(i)(ri) 



N-1 



a) cos(aJ(ri)) 



• (17) 



Equation (pTj) equals (up to an offset) the outcome of 
the classical grating. As is well-known from the grating 
equation (and as Eq. ([T7j) explicitly shows) two principal 



FIG. 2: (Color online) Left (a) - (c): density plots of 
G*^^-* (ri, r2) for two atoms versus ^(ri) and ^(r2); left (d): 
density plot of G*^'^-* (ri, r2, rs, r4) for four atoms versus ^(ri) 
and ^(r2), with ^(rs) = 7r/2 and ^(r4) = — 7r/2. Right: cuts 
through the density plots along the indicated lines, i.e., for (a) 
^(r2) = 6{ti), (b) ^(r2) = const and (c), (d) ^(r2) = -^(ri). 



maxima appear in the far-field diffraction pattern only 
if the interatomic distance is greater or equal to A (see 
Fig. [3]). By contrast, the use of the A^th order correlation 
function with the N detectors positioned according to 
Eq. ([To]) (or Eq. (p!2|) ) allows to resolve an atom- atom 
separation as small as A /A/" (or A/(A^ + 1)) (see Fig. [3]). In 
this way the A^th order correlation function G^^\ri) can 
be used to resolve and image trapped atoms separated by 
a distance d = X/N. 

Finally, we address the technical feasibility of our 
scheme. For the ability to localize atoms and adequately 
resolve optical path differences on a scale smaller than 
A we refer to [3l|, [H, [1^ . A detector of a given size s 
positioned at a distance L = |r^| in the far-field region 
(see Fig.[T]) gives rise to an angular resolution = s/L, 
i.e. to a phase resolution AS = kd cos AO. To resolve 
the modulation of the A'-th order correlation function 
G^^^(ri), a sufficient requirement is that NAS <C 27r, 
i.e. <C X/{Nd)^ which yields the condition 



Nd 

L > 5—-. 
A 



(18) 



For given A^ and d we can thus find for any detector size s 
a distance L to achieve the necessary resolution. Hereby, 
choosing the smallest L compatible with Eq. ([18]) is fa- 
vorable in order to maximize the A" photon detection 
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G(l)(ri) 



27r ^(ri) 




27r ^(ri) 



FIG. 3: (Color online) G^^\ri) and G'^^^(ri) as a function 
of ^(ri) for a chain of 4 atoms. The interval [— 27r, 27r] cor- 
responds to the maximal range of variation of ^(ri) for an 
interatomic distance d = X [see Eq. (|6])]. The dashed lines 
indicate the corresponding range in case of d = A/4. 



probability; the exact longitudinal positions of the de- 
tectors are thereby not important. In case of a gaussian 
distribution of the phases ^(r^) (i = 2, ...,A^) around 
their ideal values given by Eqs. (p!Q|) or (p!2|) with a stan- 
dard deviation a the contrast of the G*^^^ (ri)-function 
is reduced and given by e~^^ For N = 2 and 
A/" = 4, this means that a contrast of higher than 50% 
can be maintained as long as a is less than 0.8 and 
1.2, respectively. Using the set of reasonable parame- 
ters d = 5fim, Ad = 0.1 /im, 6>(ri) 30°, A(9(ri) 0.1°, 
k = 27r/800nm, Ak < 10'^ k we obtain a ^ 0.7. 

In conclusion we have shown that N photons of wave- 
length A spontaneously emitted by N atoms and coin- 
cidentally recorded by N detectors at particular posi- 
tions exhibit correlations and interference properties sim- 
ilar to classical coherent light of wavelength X/N. The 
method requires neither initially entangled states nor 
multi-photon absorption, only common single-photon de- 
tectors. 
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